livio montanaro 18/08/2022
calcolo completo direzione primaria 
calcolo del promissore al polo del significato
metodo Argoli

Il risultato è quello che vedete 
sul software morinus con l'opzione
	placidus (underthepole) 
	zodiacal usa latitudine del significatore
	aspetti del promissore al significatore.
	pds in chart:	celestial,from the pianets.

ob 	= obliquità 23.438130
lat	= latitudine luogo

dati provenienti direttamente dal vostro software

a       = significatore (0 = sole 1 = luna da tabella effemeridi)
aa      = promissore

lo[0][14] 	=longitudine ascendente
lo[0][a]	=longitudine pianeta significatore
dp      	=declinazione pianeta
arp[a]		=ascensione retta significatore
		se l'ascensione retta non è presente occorre per il significatore usare routine @cra inserendo la latitudine (dp) del pianeta,
		per calcolare ascensione retta
		con la massima precisione possibile

dh      = distanza oraria
DA      = differenza ascensionale
dr      = distanza retta
Sad/San = semiarco diurno/notturno
armc    = armc
aric    = armc + 180
qd      = quadrante significatore

ch      = circolo orario o polo significatore
ao      = ascensione obliqua
aoch    = ascensione obliqua sotto polo
dach    = DA sottopolo 


aaa     = aspectangle[aaa] = aaa = tutti gli aspetti es 0,60,90,120 ecc

comandi extra planetdance:
distance=command pd corrisponde all'angolo minore fra due punti
normalize= command pd riduce entro 360°

:psig;`calcolo significatore parte 1
u=1;
lob = lo[0][a];los=lob;`long significatore
dp = decp[a];	`declinazione significatore
lp =latp[a];	`latitudine significatore
ar = arp[a];	`a.r. pianeta significatore 
		se non presente usare routine cra con latitudine,
		per calcolare ascensione retta
		con la massima precisione possibile
DA = tan(lat)*tan(dp);
DA = asin(DA);`differenza ascensionale

Sad[u] = 90 + DA;
San[u] = 90 - DA;
dm_mc[u] = normalize(ar-armc);if (dm_mc[u] > 180) dm_mc[u] -= 360;
dm_ic[u] = normalize(ar-aric);if (dm_ic[u] > 180) dm_ic[u] -= 360;

if (normalize(-lo[0][14] + lob) > 180)
   {dh = abs(dm_mc[u])/(Sad[u]/6);} detemina se il significatore si trova sopra orizzonte
else
   {dh = abs(dm_ic[u])/(San[u]/6);} detemina se il significatore si trova sotto orizzonte
if (normalize(-lo[0][14] + lob) > 180)
   dr  = distance(armc,ar);`distanza retta
else
   dr  = distance(aric,ar);
qd = 4;
`if (dr = normalize(aric-ar)) qd = 4;
if (dr = normalize(ar-aric)) qd = 3;
if (dr = normalize(armc-ar)) qd = 2;
if (dr = normalize(ar-armc)) qd = 1;

`if (qd = 4)aoch = aric-(15*dh);` 1-3
if (qd = 3) aoch = aric+(15*dh);` 4-6
if (qd = 2) aoch = armc-(15*dh);` 7-9
if (qd = 1) aoch = armc+(15*dh);` x-12
aoch=normalize(aoch); ascensione obliqua significatore sotto il suo polo

`calcolo polo del significatore o circolo orario
DA = tan(lat)*tan(dp);
DA = asin(DA);
ch = sin(1/6*DA*dh)*1/tan(dp);
ch = atan(ch);`polo del significatore
qds = qd;`quadro dove si trova significatore (Argoli)
return;

:ppro;`promissore sotto polo del significatore parte 2
ora si calcola l'angolo del promissore alla congiunzione o agli aspetti
partendo dalla longitudine del promissore,aggiungendo o sottraendo l'angolo voluto

lob = normalize(lo[0][aa]+/-aspectangle[aaa]);`nuova longitudine del promissore

'NOTA la direzione è sempre calcolata con latitudine 0 (lp = 0) proiettata sull' equatore

lop=lob;
lp=0;@cra;calcola nuova ascensione retta sempre a latitudine 0.

dach = tan(dp)*tan(ch);	calcolo differenza ascensionale sottopolo
dach = asin(dach);
dach = abs(dach);`fondamentale
if (lat<0) dp = dp*-1;`cambio segno x temi australi

`trova direzione AO 
if (dp>0 )
   {ao = ar-dach;
    do = ar+dach;}
else
   {ao = ar+dach;
    do = ar-dach;}

if (qds = 1) delta = distance(aoch,ao);` x-12
if (qds = 2) delta = distance(aoch,do);` 7-9
if (qds = 3) delta = distance(aoch,do);` 4-6
if (qds = 4) delta = distance(aoch,ao);` 1-3

delta = angolo direzione finale

si consiglia di usare 1 grado = 1 anno

return;